NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 





(wA5A-Ca-1o4v)4 1) Ui.SCttiit'ANCici dETJhtii* 
trtPlHICAi. AND rUEOHLTICAL flODLLa UF 1'ttL 
FLkuLHG SOLkli CdHOdOSPUbHE Ai»D Tdtlri 
t'oSS I bLE itESOLUTiON (Stautoru Ult.LV.) ^7 


N ri 1 - 1 9 y o 4 


IJltO 1 <1S 




tic AjJ/MF aoi 


CSCL J 3d Gi/JZ 4 I o 4 1 






t 



DISCREPANCIES BETWEEN THEORETICAL 
AND EMPIRICAL MODELS OF THE FLARING 
SOUR CHROMOSPHERE AND THEIR 
POSSIBLE RESOLUTION 


A. Cordon Emslie 
Institute for Plasma Research 
Stanford University 


John C. Brown 
Department of Astronomy 
University of Glasgow 


<08% 


%&.%> w% 


Marcos E. Machado 
CNIE 

Buenos Aires, Argentina 


</V f* lb 


w 


a 

•> m 




National Aeronautics and Space Administration 
Grant NGL 05-020-272 


Office of Naval Research 
Contract N00014-75-C-0673 


SUIPR Report No. 820 
August 1980 


INSTITUTE FOR PLASMA RESEARCH 
STANFORD UNIVERSITY, STANFORD, CALIFORNIA 


I 


DISCREPANCIES BETWEEN THEORETICAL AND EMPIRICAL MODELS 
OF THE FLARING SOLAR CHROMOSPHERE AND THEIR POSSIBLE RESOLUTION 


A. Gordon Emslie 
Institute for Plasma Research 
Stanford University 


John C . Brown 
Department of Astronomy 
University of Glasgow 


Marcos E. Machado 
CNIE 

Buenos Aires, Argentina 


National Aeronautics and Space Administration 
Grant NGL 05-020-272 

Office of Naval Research 
Contract N00014-75-C-0673 


SUIPR Report No. 820 
August 1980 


Institute for Plasma Research 
Stanford University 
Stanford, California 


ABSTRACT 


Possible sources of tht narked dicrspancy between empirical and 
Chaoratical (paticularly a lac cr on-haa tad ) models of the solar chromosphere 
during flares are discussed. A major source of uncertainty in empirical 
models is the inhomogeneity of the spectral data on which they are baaed. 
Theoretical models Involve several possible sources of error of which 
probably the most important is neglect of the radiative coupling of upper 
and lower chromospheric regions. In addition, both types of model are 
limited by their assumption of plane-parallel, and frequently static, 
structures. 

The variety of deconvolutions involved in the construction and inter- 
comparison of empirical and theoretical models is another source of diffi- 
culty in model testing. A new procedure for studying flare energy input 
is suggested whereby the required input is derived from the empirical 
model chromosphere; this procedure is applied to the electron-heated 
case, and it is found that the integral equation defining the flare energy 
deposition rate (erg g * s can be analytically inverted to yield the 
injected electron flux energy spectrum from knowledge of the energy balance 
in the empirical atmosphere. We analyze recent empirical model results in 
this manner and compare the calculated injected electron flux spectrum with 
that needed for hard X-ray bursts in moderately large flares. This technique 
betters the usual procedure of predicting a model chromosphere from the 
hard X-ray electron flux in that it yields a measure of significance and 
uniqueness of the heating model parameters. 

Subject Headings : hydromagnetics - Sun: chromosphere - Sun* flares - 

Sun: X-rays 
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I INTRODUCTION 

A considerable effort in flare physics has recently gone into the study 
of the secondary response of the solar atmosphere to the primary flare energy 
release (see,e.g. , Canfield et al. 1980 for a review) , Essentially this field 
comprises the comparison of flare model atmospheres constructed in two dis- 
tinct ways, vis. theoretically and empirically. A theoretical model is an 
atmospheric structure predicted by solution of the equations (e.g. transport 
of energy, momentum and radiation) describing the response of the preflare 
solar plasma to a hypothetical flare energy input. An empirical model, on 
the other hand, has a structure derived from observations of the flaring 
atmosphere, by trial and error fitting of the computed model spectrum and 
subsequent model adjustment. The theoretical model may then be compared with 
the empirical or (which is not entirely equivalent - see SHI) be used to 
predict a spectrum for direct comparison with observations. Recent work 
(Machado and Emslie 1979; Lites and Cook 1979; Machado et al. 1980”“ hereafter 
MAVN) on empirical models has established marked discrepancies between these and 
theoretical models, particularly those inwhich the hypothetical energy input 
is particle beam collisions (e.g. Brown 1973, Brown, Canfield and Robertson 
1978— hereafter BCR; La3onte 1978; Machado, Emslie, and Brown 1978). Not only 
is there a quantitative discrepancy in the parameters required to describe the 
theoretical energy input (e.g. electron flux) and the observations (e.g. 

Ha brightness) of part of a particular flare (see BCR) but also qualita- 
tive differences between the general form of empirical and theoretical 
structures for any input parameters. In this paper we discuss possible 
sources of these discrepancies with a view to clarifying the most important 
areas for future work on this problem. Our dlacuasion is concerned mainly 
with models of the chromosphere, and specifically electron-heated models. 


Figure 
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The problem of the optically thin corona and tranaition region haa been 
diecusead recantly by Hachado and Ensile (1979) and that of the tempera- 
ture minimum region by Machado, Ensile, and Brown (1978). However, many 
of our remarks on the sources of discrepancy in the chromosphere are 
equally applicable to these regions and to models heated by other forms 
of energy Input, such as soft X-ray irradiation (Somov 1975; Machado 1978). 

In Figure 1 (from MAVN ) we compare the empirical models 

derived from Skylab ATM data with the electron heated (BCR) theoretical models. 
The main discrepancies are: 

(i) The BCR models dy/e higher temperatures In the mid-chromosphere 

(T * 10 4 K, h - 1100 km) and much more material with 10 4 < T < 2 * 10 4 K 

(ii) The transition region in the BCR models occurs considerably higher 
in the atmosphere than in the empirical models. Indeed for small 
electron fluxes the BCR transition region occurs higher in the 
flaring atmosphere than in the quiet sun (Machado and Emslie 1979). 

In the region h - 800-1100 km (T - 8 * 10 3 - 10 4 K) the BCR and empirical 
models have similar temperatures; to that extent the BCR models are locally 
plausible, reproducing observed Ha profiles when a moderate amount of 
horizontal inhomogeneity (see §IIc), consistent with the filling factor of the 
Ha spectrograph slit, is included. (Note however that the electron flux 
needed to obtain this agreement is ^ 1/10 of that required for thick target 
interpretation of hard X-rays from the same flare - BCR.) Nevertheless such 
a spatially localized regime as that of Ha formation clearly cannot be 
expected to provide an adequate empirical test of a model energy input. 
Furthermore (see §IIb) a localized treatment is not satisfactory theor- 
etically either. It is therefore essential to take a global view in 
modeling of the chromosphere, and indeed of the whole atmosphere. 

In ?II we consider the possible sources of the discrepancies outlined 


above and discuss the global radiative transfer problem more fully. In 
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S1II we discuss the nature of the intercomparison proceee Itself 
and how It renders estimation of uncertainties difficult. We then 
suggest a new intercomparison technique which allows estimation of these 
uncertainties; this technique involves computing the required energy 
input for a given empirical atmosphere structure, rather than constructing 
a theoretical model atmosphere baaed on a prescribed hypothetical flare 
energy input. In S1V we apply this technique to MAVN empirical flare 
data and discuss the results obtained. 
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II* SOURCES OF DISCREPANCY 

*) General 

There arc several features la common to both empirical and theoretical 
models which reader them generally unsatisfactory. First , in order to make 
the analyses tractable, most chromospheric modeling has thus far been based 
on the assumption of a steady state plane-parallel structure In hydrostatic 
equilibrium. These assumptions ell* .ate time-dependence from the equations 
and permit representation of the models purely In terms of their temperature 
versus height structure T(h) as used in Figure 1, together with a 
scaling factor defining pressure or density at a single point. However, 
they result in omission of energy transport and spectroscopic effects 
due to mass motion, time dependent conduction (see Craig and McClymont 
1976 ) and radiative transport (see Dubov 1963 ), and inhomogeneity. For 
example, if large velocity gradients were present in the atmosphere, 
then it is possible that they would radically alter the interpretation 
of spectral diagnostics in terms of a T(h) structure. (There is, however, 
little evidence for such velocity gradients in the layers of the flaring 
atmosphere under consideration [see Lites and Cook 1979 ].) Second, 
chromospheric modeling of both types involves the complex numerical Si’ution 
of the relevant radiative transfer equations, the techniques of which are 
still under development (see Vernazza, Avrett, and Loeser [VAL] 1973 , 

1976, 1980) and thus require cautious treatment, 
b) Deficiencies in Theoretical Models 

A number of factors in the electron beam energy input used by BCR warrant 
further study. First, the input is purely collisional (see Brown 1973; 

Emslie 1978) taking no account of the effects of plasma wave generation 
^rown and Melrose 1977; Vlahos and Papadopoulos 1979) nor of reverse 


currant ohaic losses (s.g. Basil* 1980) on th* height distribution of 
besa energy deposition. These ere both most important highest in the 
atmosphere end will tend to steepen the BCR T(h) profiles toward the 
satirical forma, both directly and by enhancing the role of conduction 
(see below) . Second, the effect on beam dynamics of converging magnetic 
field lines is neglected; its inclusion would enhance th* bees deposition 
at higher altitudes. Third, the BCR models were computed only for electron 
spectral index fi • 4 . Us* of the higher 5's typical of most small 
flares will also tend to steepen T(h) by depositing energy preferentially 
higher in the atmosphere (see Emslle 1978) . 

The BCR models included no thermal conduction on the grounds that 
conductive heating (or cooling) is negligible at the low temperatures 
and temperature gradients where Ha forms (see Brown 1974, 1977). We 
have checked this assumption directly for both the theoretical and empiri- 
cal models and find (see fill) that for T < 12000 K the conductive term 


is totally negligible compared to the radiative losses there so that an 
external input is needed. This is contrary to claims by Svestka (1970), 
Shmeleva and Syrovatskii (1973), and Davis, Kepple, and Strickland (1977) 
that local conductive deposition can balance chromospheric radiative 
losses during the post-flash phase of the flare. At first sight this 
also appears contrary to the statement by Machado and Emslie (1979) that 
"thermal conduction alone is sufficient to produce the observed flare 
enhancements at all atmospheric levels, in particular in the chromosphere". 
However, it turns out that in the proper, global, approach to the theoreti- 
cal modeling, it is essential to include conduction in the hotter layers 



(upper chromosphere) in order for the correct structure in the lover (Ho) 
chromosphere to be obtained. The reason is that these two regions are 
radiatively coupled, particularly by La which constitutes an Important 
energy loss frjm the upper chromosphere but a significant input into the 
lover chromosphere — Machado and Emalie (1979); MAVN. (The possibility 
of radiative heating vaa noted by But but not incorporated in their 
models.) The importance of this La "backwerming" stems from the fact 
that the energy in the La radiation comes from dissipation of thermal 
conductive flux, where the steep transition zone temperature gradient is 
suddenly reduced. In short, the Machado and Emslie (1979) statement refers 
to the fact that a thermal conductive flux in the transition cone can, 
through La backwarming , supply the energy for a significant fraction 
of the radiative losses in the upper chromosphere (i.e. down to about the 
region where Ha forms). 

The excessively high BCR transition region is attributable substantially 

to two factors. First BCR take the quiet sun model of VAL 1973 as their 

preflare state in establishing the density scale and the ambient energy input 

in their models (see Brown 1973). Use of a plage preflare model would have 

resulted in deeper transition regions, particularly in low flux models. The 

residual discrepancy, and in particular the fact that the BCR transition zone 
for low electron fluxes lies above its quiet sun level, is a numerical 

artifact. Since BCR were attempting to model the Ha formation region, and 

since they did not recognize the importance of the radiative coupling of 

upper and lower chromospheric regions, they used only a very broad grid 

representation of the upper chromosphere. In fact the grid step used 

there was of the order of the discr mcy in transition rone height noted 


by Machado ard Emslie (1979). 


c) Deficiencies In Empirical Models 


Aside from the points raised in file, end the problems of mode. unique- 
ness discussed below (fill), the mein deficiencies in empirical models lie in 
the inhomogeneity of the data used. This takes several forms. First, the 
evident inhomo gsnelty of the horltontal structure of the flare within the 
instrumental field means that the empirical models represent soma horlsontelly 
averaged flare structure. It is therefore not particularly meaningful to 
coopers these with theoretical modela representing the flare-heated areas 
(as MAVN do in Figure 1) . Rather we should compare them with 

a theoretical model diluted by the ambient atmosphere with an appropriate 
filling factor (which itself generally depends on height h according to 
the variation of flare area and of Instrumental field with h ) . The 

latter is the procedure followed by BCR in their Ha analysis. 

Second, even the best observational efforts to date ( MAVN ) 
have resulted in spectra based on very inconsistent coversge in both 

space and time. In some cases the data used in different spectral bands 
were not even from a single flare, and at best were from different parts of 
a single flare at different times. Furthermore ouch of the ATM flare cover- 
age was in the post-flash phase, so putting in doubt the relevance of the 
empirical models in testing theoretical models heated by electron beams 
which are early-stage phenomena. 

Third, the flare set analyzed by MAVN comprised five 

subflares and one Class 1, while some of the theoretical model analyses 
(e.g. BCR) have been aimed at very large events. 

III. CONVOLUTIONS AND STABILITY IN MODEL BUILDING 

An unsatisfactory aspect of current empirical models is that they give no 
error bars on T(h) so that any comment on their compatibility with theore- 


10 - 


tlcal no da Is la only qualitative. Tha quaatlon wa raally want to aak lit what 
manifold of empirical modal atmoapharaa la capable of fitting the spectral 
data within their error bars, and what theoretical models fall within this 
manifold? The first part of this question is hard to answer due to the very 
large computational effort involved In iteration to obtain even Me acceptable 
empirical model. Furthermore error bars on T(h) so obtained are almost 
certain to be rendered large and hard to interpret meaningfully because of 
the Inhomogeneity of data used. Also, in comparing theoretical end empirical 
models quantitatively, we should develop an objective criterion of acceptabi- 
lity, presumably based on some statistic such as the root mean square residual 
on the energy balance parameters (rather than simply on temperature) relative 
to the probable error in the observed (empirical) values. 

Even given a good set of homogeneous data with error bars, there are 
still fundamental limitations on the extent to which empirical models can 

test theoretical models because of the number of deconvolution processes 
involved (see Craig and Brown 1976). Specifically, in testing electron- 
heated chromospheric models, the usual procedure is (a) derive an electron 
spectrum from a noisy observed hard X-ray speetrum by deconvolution (Brown 
1971; Craig 1979); (b) predict a theoretical model atmosphere heated by 

this electron input by solving the integro-differenr lal equation of radiative 
transport; (c) derive an empirical model atmosphere structure by decon- 
volution of the integro-differential equation governing the formation of 
the observed spectrum; (d) compare (c) with (b) . Given the problems 
of error magnification in deconvolving from and of nonuniqueness in 
model-fitting to, integral equations (Craig and Brown 1976 and references 
therein), we should clearly be circumspect in rejecting or accepting 
theoretical model chromospheres heated by electrons, since their computa- 
tion depends intrinsically on the solution of three integral equations. 
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Here m suggest ths following si tarns tivs to the usual sequence above , 
which enables us to quantify uncertainties in derived quantities; 
the usual procass of forward fitting of aodels to data doea not permit 
us to do this. 

Instead of predicting an optlcal/UV spectrum from a theoretical 
model flare chromosphere with the energy input function Inferred from 
independent observations and comparing this prediction with optical and 
UV observations, we suggest reversing the procedure. That Is, we derive 
the energy input function which is consistent with the heating rates 
derived from the optlcal/UV spectral data and compare this with other 
Independent estimates of this input function. (See Machado, Ensile, 
and Brown 1978 for similar discussion of heating at the temperature minimum.) 
Specifically we propose the following procedure: 

(a) Obtain empirically that atmoi c eric temperature structure T(h) 
which best fits the spectral data, or better that range of T(h) 
which fit the data within its errors. 

(b) Evaluate the net radiative energy loss from this (these) empirical 

model (s) (see VAL 1980) and so deduce the total specific energy 
input I Q (h) (erg g^s'S required for a steady state. (I o (hl 
comprises a gross radiative loss lees a radiative input 

I^h] from backwaraing.) 

(c) Subtract from I (h) the specific conductive inpu: I - | — 4r (k [T] 4r)l 

o c L. 0 dn ah J 

(p ■ gas density) calculated from the T(h) found in (a), and also an 
estimate of the ambient energy input I a (h). The latter requires some 
assumption of how ,a‘>) changes between preflare and f .are condi- 
tions — e.g. constant 1^ per gram 'see Brown 1973). 
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(d) This luvii the energy input I f (h) to be supplied by the flaw 
riliut processes, which can than ba compared with theoretical 
heating models. The acceptability of the fit as the various para- 
meters of the theoretical model are varied gives an indication of 
the uncertainties in the values of these parameters, quantities which 
are difficult to obtain in the conventional comparison procedure 
(see BCR). It should also be noted that this procedure should 
enable a direct comparison of 1^ with I c * and 1^ at each 
height h , Indicating unambiguously the relative role of each 
process (as a function of height), which up to now has remained 
controversial. 

Ue illustrate the technique by considering chromospheric models heated by 
colliaional degradation of a beam of non-thermal electrons, injected vertically 


along uniform field lines. Expressing as a function of the Lagrangian column 
depth variable N - J" n(h')dh' rather than h , then the energy input 


arising from an injected number flux spectrum F o (E o ) (electrons 

cm -2 s _1 per unit E q ) is (Brown 1973; Eos lie 1978) 


I f (N) 


/v F (E )dE„ 

JL / 00 

®H J /3KN E od-3 KN/Ej)^ ' 


( 1 ) 


where K ■ 2 ire 4 A in the usual notation and m^ is the mass of a hydrogen 

K F o (/x) 

atom. Setting x * # v * 3KN, G(x) • *= ,-rr- and J(v) ■ I-(N) then 

o 2m* x ./3 f 


J(v) 


/’ 


G(x) dx 

( X— V ) 


( 2 ) 
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which is an Abel integral equation with solution (assuming G(x ± •) » 0 for 

all reasonable F (E )) 
o o 


G(x) 



J(v) + vJ' 



(v-x) 


i/3 




(3) 


( Courant and Hilbert 1953, pp, 158-159) • 
In terns of the original variables. 



which determines the form of the injected electron flux 

spectrum F o (E q ) from knowledge of the behavior of 1^ with N, obtained via 
steps (a) to (c) above. In the next Section we shall apply this method of 
analysis to empirical flare data from MAVN. 


Figure 2 
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IV APPLICATION TO OBSERVATIONS 

In Figure 2 we show the behevlor of I , I I end I obteined by 

o tot t c 

the procedure of VAL (1980) for model MAVN2 (see Figure 1). We use this 
model (es opposed to model MAVN1) because it corresponds to bright flare 
regions and so is more comparable with an atmosphere undergoing excitation 
by non- thermal particles, though it is still small compared to a BCR-sized 
flare. One immediately notes the following points from the Figure; 

(i) the local condutive input I fi is less than I Q by at least nine 
orders of magnitude and so is totally negligible (see remarks in 
§ 11 ). 

(ii) the ratio 1 /I „ (the contribution of backwarming to the 

r tot 

radiative energy budget) approaches unity (to within a factor 
of 2 2; see remarks in §IIb) at around N - 3*10 21 cm 2 , but 
falls off rapidly with depth thereafter (so that I * I q for 
depths greater than this). 

To derive I^(N) we further need to apply process (c), i.e. subtract 

the ambient (preflare) input I (N) and the local conductive input I (N) • 

a c 

As noted in (i) above, the latter of these is negligible. I can be 

obtained by applying the VAL (1980) method of calculating net radiative 

losses to a plage model atmosphere, preferably an empirical one. The only 

available empirical plage model is that due to Basri et al. (1979); however, 

this model suffers from serious errors in the upper chromospheric region, 

due to the assumption of local thermodynamic equilibrium for 

the Carbon atom (E. H. Avrett, private communication). Since this region 

is precisely the region of greatest interest to our present study, we do 

not consider it appropriate to use the Basri et al. (1979) model to obtain 

I (N) . We therefore have no alternative but to neglect I (N) and res- 
a a a 

trlct our attention to regions of the atmosphere in which I may be 


Figure 3 
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reasonablv taken to be negligible compered to 1^ ; hopefully future work on 
empirical plage modal r.cmosphares nay enable us to remove this restrictive 
condition. Figure 2 strongly suggests that I becomes important around 
N * 10 2 * cm* 2 (note the upturn in the graph of 1 (H) at this point); this 

V 

belief is somewhat strengthened by application of the VAL (1980) radiative 
loss calculation to the Basri et al. (1979) model, which shows that radiative 
losses due to Mg and Ca lines start to become very important at just this 
depth. 

In Figure 3 we show the result of applying equation (4) to the I q 

(= I. , since we are neglecting I ) values of Figure 2. To avoid the 
C s 

singularity at the lover integral limit in our computations we first trans- 
formed (4) to the equivalent expression 


F (E ) 
o o 



dx 


(5) 


(see Brown 1971). In the evaluation of this integral there is a truncation 

error at the upper limit, since we only have data up to - OKN^^/E 2 -!) 2 / 3 , 

where is the largest column density for which data is available. This 

truncation error will only seriously affect the determination of F o (E q ) near 

E = (3KN )H , which is well outside the range of E values for which a 
o max ’ o 

spectrum is to be determined (see below); thus it does not affect the results 

to follow. Note the hardening in the electron spectrum above around 300 keV, 

22 -2 

corresponding to N = 10 cm — cf . the lower integral limit in equation (4) . 

This is attributable in part to our neglect of I , as noted above, but is 

A 

certainly also in part representative cf the strong heating observed near the 
temperature minimum which cannot be explained by non-thermal electron bombard- 
ment (Machado, Emslie, and Brown 1978). Shown in the Figure is the best-fit 

power-law F (E ) ■ F E ^ to the (E ,F ) points, where only E values to the 
r oo ooo o’ o r o 

left of the spectral bend have been used (see remarks above on the limits of 


applicability of our method). This best-tit line has the form 
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F (B ) • 0.53) 

O 0 


e -( 6.21 ±0.24) 
o 


electrons cm” Z 8 _l k«v” 1 , 


( 6 ) 


with a regression coefficient r ■ 0.997 . 

For comparison with hard X-ray burst anal y a la it la more convenient to 
axpraaa In terms of the total Injected energy flux above 20 keV (any), 

via. 

F J# (erg cm"V l ) - 1.6 * I0~'f T q (Z q ) E q dE Q . (7) 


Using equation (6) in equation (7) ve obtain, for the Injected electron flux 
parameters: 


log F - 12.01 ± 0.53 
20 

6 - 6.21 ± 0.24 


( 8 ) 


From these results and their uncertainties ve conclude that heating 
by an electron beam with a power law injection spectrum fits the observa- 
tions extremely well over the atmospheric domain concerned; however, the 
scatter in the few data points used prevents us from deriving the required 
electron flux F 2Q with an accuracy better than about half an order of 
magnitude. Ue also note that the value of 6 in eq. (8) is much larger 
than the 5-4 used by BCR. This explains the higher temperatures (compared 
to the MAVN models) found by BCR in the Ha formation region (II); see 
comments in ill on the effect of different 5 on the variation of the 
energy deposition rate with depth. With this in mind, we note, however, 
that the electron injection parameters (8) are entirely consistent with 

those required for hard X-ray bursts in large events (Hoyng, Brown, and 

—I 

van Beek 1976), viz. F,,, 2 3 * 10* 1 erg cm s (for plausible injection 
areas A * 10 18 - 10 1 - era-) and 5=4-6. This is encouraging, although 
it must be noted that the above illustrative analysis is based on a very 
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Halted range of points, vis. a factor of ^5 in the column density N 

and so a factor of ^2 in the electron energy E q . In addition, the 

error on F 2q quoted in eq. (8) is based on the residuals of the 

(E ,F ) points about the best fit line (Figure 3) and might be substantially 
o o 

larger if the uncertainties in the individual points exceed these residuals; 
these uncertainties are, however, very difficult to estimate (E.H. Avrett, 
private coMunication) . It is to be hoped that more work on preflare 
structures, and the acquisition of more extensive data sets, will enable 
the present method to be used as a better flare diagnostic technique. 


Wa thank E.H. Avratt and R. Loasar for perforaing tha t ima-consuaing computational 
work involved in dariving tha nat radiative loss rates shown in Figure 2. 
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FIGORE CAPTIONS 

Flgurt It Tsmpsraturs versus height structures in empirical model chromospheres 
for smell (MAVN i) end larger (MAVN 2) flares constructed from 
Sky lab data by Machado at al^ (1980), compared with the structures 
predicted by Brown, Canfield, and Robertson (1978) for chromospheres 
heated by electron beams with power law spectral Index 4 and energy 
fluxes (above 20 keV) of 10 9 (BCR 9) and 10 10 (BCR 10) ergs cm* 2 s* 1 
respectively. The column density values at the top of the Figure 
correspond to model MAVN 2 (cf. Figure 2). 

Figure 2: Energy sources and sinks as functions of column density N 

inferred from the empirical MAVN 2 flare structure. 1 (N) is the 

o 

net radiative loss (ergs g~ l s* 1 ), comprising a gross radiative 
loss term l tQt and a (negative) backwarming term . I C (N) 
is the conductive deposition inferred from the temperature 
gradients in the atmosphere. Note the importance of backwarming 
effects for N £ 10 22 cm* 2 and the negligible amount of local 
conductive heating. 

Figure 3: Injected electron flux spectrum obtained by inversion of the I 0 (N) 

values in Figure 2 (see equation [4]). Beyond E - 300 keV 

22 -2 

(corresponding — see equation [4]— to N s 10 cm ) the spectrum 
changes form dramatically; this is partly real and partly due 
to the neglect of preflare heating terms in the flare energy 
equation — see discussion in text (Sill). Below this break point, 
the spectrum is well fit by a power law with the parameters shown; 
these can be compared to independently derived electron flux 
parameters, obtained from analysis of hard X-rcy spectre (see 
equation [8]). 
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